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ABSTRACT 

HR  8799  is  an  hFO  mA5  y  Doradus-,  X  Bootis-,  Vega-type  star  best  known  for  hosting  four  directly  imaged  candidate 
planetary  companions.  Using  the  CHARA  Array  interferometer,  we  measure  HR  8799’s  limb-darkened  angular 
diameter  to  be  0.342  ±  0.008  mas  (an  error  of  only  2%).  By  combining  our  measurement  with  the  star’s  parallax 
and  photometry  from  the  literature,  we  greatly  improve  upon  previous  estimates  of  its  fundamental  parameters, 
including  stellar  radius  (1.44  ±  0.06  R0),  effective  temperature  (7193  ±  87  K,  consistent  with  F0),  luminosity 
(5.05  ±  0.29  Lq),  and  the  extent  of  the  habitable  zone  (HZ;  1.62-3.32  AU).  These  improved  stellar  properties  permit 
much  more  precise  comparisons  with  stellar  evolutionary  models,  from  which  a  mass  and  age  can  be  determined, 
once  the  metallicity  of  the  star  is  known.  Considering  the  observational  properties  of  other  X  Bootis  stars  and 
the  indirect  evidence  for  youth  of  HR  8799,  we  argue  that  the  internal  abundance,  and  what  we  refer  to  as  the 
effective  abundance,  is  most  likely  near  solar.  Finally,  using  the  Yonsei-Yale  evolutionary  models  with  uniformly 
scaled  solar-like  abundances,  we  estimate  HR  8799’s  mass  and  age  considering  two  possibilities:  1.5167oo34  M0 
and  3347132  Myr  if  the  star  is  contracting  toward  the  zero-age  main  sequence  or  1.5134oo24  and  9043580I  Myr 
if  it  is  expanding  from  it.  This  improved  estimate  of  HR  8799’s  age  with  realistic  uncertainties  provides  the  best 
constraints  to  date  on  the  masses  of  its  orbiting  companions,  and  strongly  suggests  they  are  indeed  planets.  They 
nevertheless  all  appear  to  orbit  well  outside  the  HZ  of  this  young  star. 
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angular  resolution  -  techniques:  interferometric 
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1.  INTRODUCTION 

HR  8799  (HD  218396,  HIP  1 14189)  is  a  y  Doradus-type  star, 
a  class  characterized  by  pulsations  on  timescales  longer  than  the 
typical  <5  Scuti  stars,  their  neighbors  on  the  Hertzsprung-Russell 
(H-R)  diagram  (Breger  1979;  Zerbi  et  al.  1999;  Kaye  et  al.  1999; 
Handler  1999).  HR  8799  is  also  a  X  Bootis-type  star,  which  is 
a  spectroscopically  defined  group  of  non-magnetic,  chemically 
peculiar,  Population  I  A-  and  early  F-type  stars  that  show  metal 
deficiencies  in  Fe-peak  elements  while  showing  solar  to  slightly 
over  solar  abundances  of  C,  N,  O,  and  S  (Venn  &  Lambert 
1990;  Paunzen  et  al.  1998;  Gray  &  Kaye  1999).  In  addition,  HR 
8799  is  a  “Vega-like”  star,  possessing  excess  infrared  emission 
longward  of  ~20  // m  (Sadakane  &  Nishida  1986;  Su  et  al.  2009), 
attributed  to  thermal  emission  from  a  debris  disk. 

Aside  from  these  three  distinguishing  characteristics,  interest 
in  HR  8799  was  amplified  in  2008  with  the  discovery  of  three 
planet-like  companions  at  projected  separations  of  24,  38,  and 
68  AU  (Marois  et  al.  2008).  From  comparisons  of  the  observed 
fluxes  to  the  predictions  of  evolutionary  models,  Marois  et  al. 
estimated  companion  masses  of  5-11  Mjupjter  for  companion  b 
and  7-13  MJupiter  for  companions  c  and  d,  each.  Because  the 
masses  of  the  companions  depend  on  the  age  of  the  system, 
Marois  et  al.  used  four  lines  of  reasoning  to  estimate  the  host 
star’s  age:  (1)  comparing  the  star’s  galactic  space  motion  to 


other  young  stars  in  the  solar  neighborhood,  (2)  placing  the  star 
on  a  color-magnitude  diagram,  (3)  characterizing  the  typical 
ages  of  X  Boo  and  y  Dor  stars,  and  (4)  the  assumption  that  HR 
8799  has  a  massive  debris  disk,  which  is  typically  only  found 
to  be  associated  with  young  stars  (<500  Myr).  Based  on  these 
characteristics,  they  constrained  an  age  range  between  30  and 
160  Myr.  A  fourth  companion  was  imaged  two  years  later  with 
an  estimated  mass  of  7432  Mjup;tel-  if  the  system  is  30  Myr  old 
and  10  ±  3  Mjupiter  if  60  Myr  old  (Marois  et  al.  2010). 

Marley  et  al.  (2012)  fit  atmospheric  and  evolution  models  to 
the  data  to  determine  the  masses,  radii,  temperatures,  gravities, 
and  cloud  properties  for  the  planets.  They  found  masses  of  26 
Tfjupiter  for  planet  b  and  8-1 1  Mjupiter  for  planets  c  and  d,  though 
they  acknowledge  that  their  model  fit  for  b’s  mass  is  not  likely  to 
be  the  true  value.  They  also  determined  an  age  range  of  360  Myr, 
40-100  My,  and  30-100  Myr  for  each  planet,  respectively. 

The  ages  adopted  by  Marois  et  al.  (2008,  2010)  have  become 
a  topic  of  some  debate.  The  star’s  age  is  of  vital  interest  because 
it  is  directly  linked  to  the  masses  of  the  companions.  If  the  star  is 
young,  its  companions  are  brighter  and  their  inferred  masses  are 
lower — i.e.,  planets — than  if  the  star  is  older,  in  which  case  its 
more  massive  companions — i.e.,  brown  dwarfs — have  cooled 
significantly  and  could  be  mistaken  for  planets  (Moro-Martln 
et  al.  2010).  For  example,  if  HR  8799  is  as  old  as  the  Hyades 
(625  Myr;  Perryman  et  al.  1998),  all  companions  would  be 


1 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

21  NOV  2012 


2.  REPORT  TYPE 


4.  TITLE  AND  SUBTITLE 

The  Chara  Array  Angular  Diameter  of  HR  8799  Favors  Planetary 
Masses  for  Its  Imaged  Companions 

6.  AUTHOR(S) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory, Remote  Sensing  Division, 4555  Overlook 
Avenue  SW, Washington, DC, 20375 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


3.  DATES  COVERED 

00-00-2012  to  00-00-2012 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROIECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

HR8799  is  an  hFO  mA5  &#947;  Doradus-,  &#955;  Bootis-, Vega-type  star  best  known  for  hosting  four 
directly  imaged  candidate  planetary  companions.  Using  the  CHARA  Array  interferometer,  we  measure 
HR  8799?s  limb-darkened  angular  diameter  to  be  0.342  ?  0.008  mas  (an  error  of  only  2%).  By  combining 
our  measurement  with  the  star?s  parallax  and  photometry  from  the  literature,  we  greatly  improve  upon 
previous  estimates  of  its  fundamental  parameters  including  stellar  radius  (1.44  ?  0.06  R ),  effective 
temperature  (7193  ?  87  K,  consistent  with  F0),  luminosity  (5.0570.29  L  ),  and  the  extent  of  the  habitable 
zone  (HZ;  1.6273.32  AU).  These  improved  stellar  properties  permit  much  more  precise  comparisons  with 
stellar  evolutionary  models,  from  which  a  mass  and  age  can  be  determined  once  the  metallicity  of  the  star 
is  known.  Considering  the  observational  properties  of  other  &#955;  Bootis  stars  and  the  indirect  evidence 
for  youth  of  HR  8799,  we  argue  that  the  internal  abundance,  and  what  we  refer  to  as  the  effective 
abundance,  is  most  likely  near  solar.  Finally,  using  the  Yonsei-Yale  evolutionary  models  with  uniformly 
scaled  solar-like  abundances,  we  estimate  HR  87997s  mass  and  age  considering  two  possibilities: 
1.516+0.038  &#8722;0.024  M  and  33+7  &#8722;13.2  Myr  if  the  star  is  contracting  toward  the  zero-age 
main  sequence  or  1.513+0.023  &#8722;0.024  M  and  90+381  &#8722;50  Myr  if  it  is  expanding  from  it.  This 
improved  estimate  of  HR  87997s  age  with  realistic  uncertainties  provides  the  best  constraints  to  date  on  the 
masses  of  its  orbiting  companions,  and  strongly  suggests  they  are  indeed  planets.  They  nevertheless  all 
appear  to  orbit  well  outside  the  HZ  of  this  young  star. 

15.  SUBIECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF 

ABSTRACT 

OF  PAGES 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

15 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


The  Astrophysical  Journal,  761:57  (15pp),  2012  December  10 


Baines  et  al. 


Table  1 

HR  8799  System  Parameters  from  the  Literature 


^fstar 

(M0) 

Age 

(Myr) 

Method  Used 

Reference 

1.47  ±0.30 

Spectral  synthesis  &  spectrophotometry 

Gray  &  Kaye  (1999) 

50-1128 

H-R  diagram  placement 

Song  et  al.  (2001) 

30 

Stellar  kinematics,  H-R  diagram  placement  vs.  isochrones 

Zuckerman  &  Song  (2004) 

20-150 

Local  Association  membership 

Moor  et  al.  (2006) 

590 

IR  excess 

Chen  et  al.  (2006) 

30 

IR  excess,  H-R  diagram  placement 

Rhee  et  al.  (2007) 

1.5 

30-160 

Multiple  methods 

Marois  et  al.  (2008) 

1 .2—1.6 

<100  Myr 

Dynamical  stability  analysis 

Gozdziewski  &  Migaszewski  (2009) 

£50 

Observational  data  analysis 

Reidemeister  et  al.  (2009) 

1.5 

<100  Myr 

Dynamical  stability  analysis 

Fabrycky  &  Murray-Clay  (2010) 

~  100  Myr 

Disk  inclination 

Moro-Mart£n  et  al.  (2010) 

1.32-1.33“ 

1123-1623 

Asteroseismology 

Moya  et  al.  (2010a) 

1.44-1.45b 

26M30 

Asteroseismology 

Moya  et  al.  (2010a) 

1.32 

1126-1486 

Asteroseismology 

Moya  et  al.  (2010a) 

30 

Columba  Association  membership 

Doyon  et  al.  (2010) 

30 

Columba  Association  membership 

Zuckerman  et  al.  (2011) 

30-100 

Direct  imaging 

Currie  et  al.  (2011) 

30-300 

Atmospheric/evolution  model  fitting 

Marley  et  al.  (2012) 

30-155 

Dynamical  stability  analysis 

Sudol  &  Haghighipour  (2012) 

Notes. 

“  Solution  1 . 
b  Solution  2. 


Table  2 

Calibrator  Star  Properties 


HD 

Spectral 

Type 

V 

(mag) 

K 

(mag) 

E(B  -  V) 

^photometric 

(mas) 

213617 

FI  V 

6.42 

5.58  ±  0.02 

0.03 

0.288  ±  0.006 

214698 

A2  V 

6.33 

6.21  ±  0.02 

0.05 

0.189  ±0.005 

219487 

F5  V 

6.60 

5.54  ±  0.02 

0.05 

0.304  ±  0.006 

Notes.  Spectral  types  are  from  SIMBAD;  V  magnitudes  are  from  Mermilliod 
(1991)  except  for  HD  214698,  which  is  from  Perryman  &  ESA  (1997).  No 
errors  were  listed  so  we  assigned  errors  of  ±0.01  mag;  K  magnitudes  are  from 
Cutri  et  al.  (2003). 


more  massive  than  13.6  Mjupjter,  and  thus  more  appropriately 
called  brown  dwarfs  according  to  popular  convention.  One 
example  of  the  effect  of  presumed  ages  on  masses  involved  the 
announcement  by  Lafreniere  et  al.  (2010)  of  a  directly  imaged 
planetary  companion  to  1RXS  J160929. 1-2 10524  in  the  Upper 
Scorpius  association,  which  relied  on  an  assumed  cluster  age 
of  5  Myr.  Soon  afterward,  Pecaut  et  al.  (2012)  determined  an 
older  age  for  the  cluster:  1  U'2  Myr.  In  this  case,  the  imaged 
companion  has  a  mass  of  ~14  Mjupiter,  putting  it  just  above  the 
planetary  mass  limit. 

Moya  et  al.  (2010a)  refute  the  Marois  et  al.  (2008)  young 
age  of  HR  8799  on  all  four  counts.  They  point  out  that 

(1)  space  motions  of  young  disk  stars  are  often  inconclusive, 

(2)  the  uncertain  internal  metallicity  of  HR  8799  makes  com¬ 
parisons  with  evolutionary  models  and  other  stellar  clusters  un¬ 
reliable,  (3)  that  A  Boo  and  y  Dor  type  stars  have  a  broad  range 
of  ages,  and  (4)  that  debris  disks  are  “highly  chaotic  and  un¬ 
predictable,”  and  cannot  be  used  to  estimate  stellar  age.  Moya 
et  al.  (2010a)  instead  determined  HR  8799’s  age  by  modeling 
y  Dor  pulsations  detected  by  Zerbi  et  al.  (1999)  and  adopting 
spectroscopic  parameters  by  Allende  Prieto  &  Lambert  (1999). 
Although  primarily  limited  by  the  unknown  inclination  of  the 


Table  3 

Calibrated  Visibilities 


Date 

Calibrator 

Spatial  Freq 

V2 

aV2 

Used 

(10s  rad-1) 

2010  Aug  25  HD  213617 

2.223 

0.861 

0.081 

2.245 

0.768 

0.069 

2.267 

0.785 

0.069 

2.289 

0.807 

0.076 

2.311 

0.804 

0.089 

2.332 

0.729 

0.094 

2.353 

0.811 

0.138 

2.375 

0.844 

0.146 

2.397 

0.924 

0.142 

2.419 

0.976 

0.135 

2.442 

0.892 

0.092 

2.464 

0.834 

0.066 

2.486 

0.788 

0.059 

2.507 

0.764 

0.057 

2.528 

0.804 

0.063 

2.550 

0.742 

0.059 

2.571 

0.728 

0.067 

2.593 

0.713 

0.084 

2.613 

0.771 

0.104 

2.634 

0.840 

0.120 

2.655 

0.784 

0.103 

2.677 

0.827 

0.096 

2.699 

0.801 

0.077 

Notes.  All  data  were  taken  using  the  S2-W2  baseline  (177  m). 

(This  table  is  available  in  its  entirety  in  a  machine-readable  form  in  the  online 
journal.  A  portion  is  shown  here  for  guidance  regarding  its  form  and  content.) 


star’s  rotation  axis,  they  find  a  preferred  age  of  between  1 . 1  and 
1 .5  Gyr;  again,  this  age  implies  that  the  companions  to  HR  8799 
are  brown  dwarfs. 

Moro-Martfn  et  al.  (2010)  subsequently  challenged  the  claim 
made  by  Moya  et  al.  (2010a),  finding  an  age  of  the  system 
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Table  4 

HR  8799  Angular  Diameter  Measurements 

Date 

Calibrator 

#UD 

#LD 

HD 

(mas) 

(mas) 

2010  Aug  25 

213617 

0.297  ±0.021 

0.309  ±0.021 

219487 

0.381  ±  0.015 

0.397  ±  0.016 

All 

0.357  ±  0.016 

0.372  ±  0.016 

2010  Sep  7 

213617 

0.342  ±  0.015 

0.356  ±  0.014 

219487 

0.375  ±  0.010 

0.391  ±  0.010 

All 

0.359  ±  0.009 

0.373  ±  0.010 

2010  Sep  8 

213617 

0.302  ±0.018 

0.315  ±  0.019 

219487 

0.313  ±  0.012 

0.326  ±  0.013 

All 

0.308  ±  0.011 

0.321  ±  0.012 

2011  Oct  20 

213617 

0.320  ±  0.013 

0.333  ±  0.014 

214698 

0.349  ±  0.017 

0.363  ±  0.018 

219487 

0.314  ±0.015 

0.327  ±  0.015 

All 

0.331  ±  0.009 

0.345  ±  0.010 

2011  Oct  21 

213617 

0.336  ±  0.010 

0.350  ±  0.010 

214698 

0.320  ±  0.012 

0.333  ±  0.013 

219487 

0.332  ±0.016 

0.346  ±  0.017 

All 

0.330  ±  0.019 

0.343  ±  0.021 

2011  Oct  22 

213617 

0.308  ±  0.010 

0.321  ±  0.011 

214698 

0.308  ±  0.014 

0.320  ±  0.015 

219487 

0.332  ±  0.008 

0.346  ±  0.009 

All 

0.328  ±  0.023 

0.341  ±  0.022 

2011  Sep  30 

214698 

0.323  ±  0.009 

0.338  ±  0.008 

All  dates 

213617 

0.327  ±  0.008 

0.341  ±  0.008 

2011  Aug  data  only 

0.326  ±  0.008 

0.340  ±  0.008 

All  dates 

214698 

0.327  ±  0.009 

0.341  ±  0.010 

All  dates 

219487 

0.343  ±  0.008 

0.358  ±  0.008 

2011  Aug  data  only 

0.328  ±  0.009 

0.342  ±  0.009 

All  dates 

All 

0.333  ±  0.007 

0.347  ±  0.007 

Final  fit:  2011  data 

All 

0.327  ±  0.007 

0.342  ±  0.008 

Notes.  Figure  10  shows 

a  graphical  view 

of  these  values. 

Data  were  obtained 

using  the  S2-W2  baseline  at  177  m  for  all  nights  except  for  201 1  Sep  30,  which 

were  obtained  using  the  Sl-El  baseline  at  331  m. 

Table  5 

HR  8799  Photometry 

Band 

Average 

No.  of 

(mag) 

Measurements 

U 

6.191  ±  0.016 

2 

B 

6.235  ±  0.016 

2 

V 

5.980  ±  0.016 

2 

U 

7.471  ±  0.025 

7 

V 

6.468  ±  0.025 

7 

b 

6.143  ±  0.027 

7 

y 

5.962  ±  0.026 

7 

j 

5.383  ±  0.027 

1 

H 

5.280  ±  0.018 

1 

K 

5.240  ±  0.018 

1 

of  ~  1 00  Myr  after  considering  the  limiting  case  for  orbital 
stability;  at  these  relatively  small  separations,  massive  com¬ 
panions  would  not  be  dynamically  stable  on  long  timescales. 
They  used  the  double  4:2:1  mean  motion  resonance  configura¬ 
tion  presented  by  other  authors  (e.g.,  Reidemeister  et  al.  2009; 
Gozdziewski  &  Migaszewski  2009;  Fabrycky  &  Murray-Clay 
2010)  as  a  way  to  ensure  dynamical  stability  by  avoiding  close 
encounters.  Moro-Martm  et  al.  used  evolutionary  models  to  pre¬ 
dict  ages  from  150  to  350  Myr  and  conclude  that  the  dynamical 
state  of  the  system  favors  an  age  of  ~  100  Myr,  which  places  the 
companions  just  shy  of  the  brown  dwarf  regime.  However,  they 


Table  6 

HR  8799  Stellar  Parameters 


Parameter 

Value 

Reference 

V  magnitude 

5.98  ±0.01 

Mermilliod  (1991) 

K  magnitude 

5.24  ±  0.02 

Cutri  et  al.  (2003) 

7 r  (mas) 

25.38  ±  0.70 

van  Leeuwen  (2007) 

Distance  (pc) 

39.40  ±  1.09 

van  Leeuwen  (2007) 

Vx 

0.49  ±  0.02 

Claret  &  Bloemen  (201 1) 

Ay 

0.00  ±  0.01 

This  work;  PED  fit 

Fbol  (10— 7  erg  s-1  cm-2) 

1.043  ±0.012 

This  work;  PED  fit 

^eff, estimated  (K) 

7211  ±90 

This  work;  PED  fit 

#UD  (mas) 

0.327  ±  0.008 

This  work 

#ld  (mas) 

0.342  ±  0.008 

This  work 

^linear  (^o) 

1.44  ±0.06 

This  work 

Teff  (K) 

7193  ±  87 

This  work 

L(L0) 

5.05  ±  0.29 

This  work 

If  the  star  is  contracting  toward  the  ZAMS 


Age  (Myr) 

33+713  2  This  work 

Mass  (M0) 

1.51 6^  024  This  work 

If  the  star  is  expanding  from  the  ZAMS 

Age  (Myr) 

9(C  5Q1  This  work 

Mass  (M0) 

1.51 3^.o  024  Thi  s  work 

point  out  that  the  average  temperature  of  7350  K  and  luminosity 
of  5  Lq  from  Sadakane  (2006)  and  Moya  et  al.  (2010b)  imply 
an  age  for  the  star  that  is  very  young  when  compared  to  other  X 
Boo  stars. 

In  addition  to  these  recent  attempts  to  estimate  HR  8799’s 
age,  many  other  efforts  have  been  made  to  determine  its  age  in 
the  last  decade.  The  resulting  values  vary  from  20  Myr  up  to 
1.5  Gyr.  In  Table  1,  we  summarize  both  the  attempts  described 
above  and  other  efforts.  A  variety  of  methods  were  utilized, 
including  inspecting  the  star’s  proper  motion,  comparing  its 
position  on  an  H-R  diagram  with  isochrones,  inspecting  the  star 
for  infrared  excess,  asteroseismology,  astrometry,  spectroscopy, 
dynamical  stability  analysis,  disk  symmetry,  direct  imaging,  and 
group  membership. 

Here,  we  seek  to  greatly  improve  our  understanding  of  HR 
8799’s  age  by  using  interferometric  observations  to  directly 
measure  its  angular  diameter.  This  value,  when  combined  with 
the  HIPPARCOS  parallax,  yields  the  star’s  physical  radius.  In 
combination  with  HR  8799’s  bolometric  flux  determined  from 
broadband  photometry,  the  stellar  luminosity  and  effective  tem¬ 
perature  are  precisely  determined.  Using  these  results,  we  pro¬ 
vide  new  age  and  mass  estimates  based  on  comparisons  with 
the  stellar  evolutionary  models,  adopting  a  metallicity  most  ap¬ 
propriate  for  X  Boo  stars  and  HR  8799  in  particular.  Specifi¬ 
cally,  Section  2  describes  the  interferometric  observations  and 
calibrator  star  selection;  Section  3  discusses  the  visibility  mea¬ 
surements  and  how  stellar  parameters  were  calculated,  including 
angular  diameter,  radius,  luminosity,  and  temperature;  Section  4 
explores  the  physical  implications  of  the  interferometric  obser¬ 
vations;  and  Section  5  summarizes  our  findings. 

2.  INTERFEROMETRIC  OBSERVATIONS 

Observations  were  obtained  using  the  Center  for  High 
Angular  Resolution  Astronomy  (CHARA)  Array,  a  six  ele¬ 
ment  optical-infrared  interferometer  located  on  Mount  Wilson, 
California  (ten  Brummelaar  et  al.  2005).  All  observations  used 
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Figure  1.  #ld  fit  for  HR  8799  on  2010  August  25.  The  diamonds  and  vertical  lines  are  the  measured  visibilities  and  their  associated  errors,  the  solid  line  is  the  best-fit 
LD  diameter,  the  dashed  line  is  the  lcr  error,  and  the  dotted  line  is  the  2 o  error. 

the  Precision  Astronomical  Visible  Observations  (PAVO)  beam 
combiner  with  a  spectral  bandpass  of  ~650-800  nm.  For  a  de¬ 
scription  of  the  PAVO  instrument  and  data  reduction  procedure, 
see  Ireland  et  al.  (2008).  We  observed  HR  8799  over  seven  nights 
spanning  two  years  with  the  S2-W2  and  Sl-El  telescope  pairs 
with  maximum  baselines  of  177  m  and  331  m,  respectively.7 

We  interleaved  calibrator  star  and  HR  8799  observations  so 
that  every  target  was  flanked  by  calibrator  observations  made 
as  close  in  time  as  possible,  which  allowed  us  to  convert  instru¬ 
mental  target  and  calibrator  visibilities  to  calibrated  visibilities 
for  the  target.  Calibrators  are  stars  with  predicted  diameters  that 
are  significantly  smaller  and  in  close  proximity  in  the  sky  to  the 
target  star.  Reliable  calibrators  were  chosen  to  be  single  stars 
with  angular  diameter  estimates  <0.3  mas  so  they  were  nearly 
unresolved  on  the  baseline  used,  which  meant  uncertainties  in 
each  calibrator’s  diameter  did  not  affect  the  target’s  diameter 
calculation  as  much  as  if  the  calibrator  star  had  a  larger  angular 
size. 


1  The  three  arms  of  the  CHARA  Array  are  denoted  by  their  cardinal 
directions:  “S,”  “E and  “W”  are  south,  east,  and  west,  respectively.  Each  arm 
bears  two  telescopes,  numbered  “1”  for  the  telescope  farthest  from  the  beam 
combining  laboratory  and  “2”  for  the  telescope  closer  to  the  lab.  The 
“baseline”  is  the  distance  between  the  telescopes. 


In  order  to  estimate  the  reddening  of  each  calibrator  star,  we 
obtained  the  Tycho  {Bj  —  Vt)  color  from  Perryman  &  ESA 
(1997)  and  converted  to  (B  —  V)  using  the  table  in  Bessell 
(2000).  We  then  compared  the  observed  (B  —  V)  value  with  the 
list  of  intrinsic  (B  —  V)  colors  as  a  function  of  spectral  type  given 
by  Schmidt- Kaler  (1982)8  to  arrive  at  an  estimate  of  E(B  —  V), 
and  adopted  the  reddening  law  described  in  O’Donnell  (1994) 
to  deredden  the  observed  magnitudes.  The  photometric  angular 
diameters  were  then  determined  using  the  relationship  between 
the  (V  —  K)  color  and  log  0LD  from  Kervella  et  al.  (2004).  The 
error  in  the  diameters  is  due  to  the  relative  calibration  error 
stated  in  the  Kervella  et  al.  (2004)  paper  as  well  as  errors  in 
photometry  measurements.  Table  2  lists  the  input  photometry 
and  resulting  photometrically  estimated  angular  diameters  of 
the  calibrator  stars. 

3.  RESULTS 

3. 1.  Angular  Diameter  Measurement 

The  observed  quantity  of  an  interferometer  is  fringe  contrast 
or  “visibility,”  more  formally  defined  as  the  correlation  of  two 


The  table  was  obtained  from  http://obswww.unige.ch/gcpd/mk01bv.html 
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Figure  2.  fit  for  HR  8799  on  2010  September  7.  The  symbols  are  the  same  as  in  Figure  1. 


wave  fronts  whose  amplitude  is  the  visibility  squared  ( V 2), 
which  is  fit  to  a  model  of  a  uniformly  illuminated  disk  (UD) 
that  represents  the  face  of  the  star.  Diameter  fits  to  V2  were 
based  upon  the  UD  approximation  given  by  V  2  =  ([2J[  (x)]/ x)1 , 
where  J\  is  the  first-order  Bessel  function  and  x  =  n  B(\-nk[ , 
where  B  is  the  projected  baseline  at  the  star’s  position,  6*ud  is  the 
apparent  UD  angular  diameter  of  the  star,  and  k  is  the  effective 
wavelength  of  the  observation  (Shao  &  Colavita  1992).  A  more 
realistic  model  of  a  star’s  disk  involves  limb-darkening  (LD), 
and  the  relationship  incorporating  the  linear  LD  coefficient  /x; 
(Hanbury  Brown  et  al.  1974)  is 


V2 


-2 


X 


(1  -  Bx) 


J\(x) 

X 


t/2  Jy 2(X) 

X 3/2 


(1) 


Table  3  lists  the  date  of  observation,  the  calibrator  used,  k/B, 
the  calibrated  visibilities  (V2),  and  errors  in  V2  ( crV 2). 

The  LD  coefficient  was  obtained  from  Claret  &  Bloemen 
(20 1 1 )  after  adopting  the  effective  temperature  (Teff)  and  surface 
gravity  (log  g )  values  from  Allende  Prieto  &  Lambert  (1999), 
which  were  7586  K  and  4.35,  respectively.  The  resulting  UD  and 


LD  angular  diameters  are  listed  in  Table  6.  Figures  1-8  show 
the  LD  diameter  fits  for  HR  8799  by  night  and  by  calibrator  and 
Figure  9  shows  all  the  data  combined. 

Table  4  lists  the  resulting  angular  diameter  fits  for  each  night 
using  each  calibrator,  and  Figure  10  shows  a  graphical  version 
of  the  results.  Though  there  is  some  scatter  in  the  diameter  fit 
from  each  night/calibrator  combination,  the  scatter  from  the 
2011  data  is  less  pronounced  than  from  2010  data,  which  is 
seen  in  the  diameter  fits  shown  in  Figures  4-7.  In  particular,  for 
observations  obtained  in  2010,  the  data  exhibit  sinusoidal-like 
variations  about  the  best  angular  diameter  fit;  the  variations  are 
not  seen  in  2011  data.  This  is  almost  certainly  due  to  coating 
asymmetries  between  CHARA  Array  telescopes  (particularly 
overcoated  silver  versus  aluminum)  that  were  present  in  2010 
but  removed  in  2011.  Attempts  were  made  to  search  for 
polarization  effects  that  would  cause  these  residuals,  but  no 
conclusive  evidence  was  found.  Visibilities  were,  however, 
never  measured  in  full  Stokes  parameters  in  a  configuration 
that  showed  these  residuals. 

Because  the  201 1  data  show  none  of  the  sinusoidal  residuals, 
our  final  0ld  incorporates  the  data  from  2011  only.  This  yielded 
an  angular  diameter  of  0.341  ±  0.008  mas.  When  the  data 
from  all  nights  and  using  all  calibrators  are  fit  together,  0LD  is 
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0.347  ±  0.007  mas,  which  is  within  the  uncertainty  of  the 
adopted  value. 

The  uncertainty  for  0LD  was  calculated  as  described  in  the 
supplemental  material  from  Derekas  et  al.  (2011):  for  each  one 
of  10,000  Monte  Carlo  simulated  Gaussian  distributions,  we 
calibrated  the  instrumental  visibilities  and  fit  an  angular  diam¬ 
eter  to  the  calibrated  data  using  a  least-squares  minimization. 
We  accounted  for  random  errors  by  adding  random  numbers 
generated  from  the  empirical  covariance  matrix  scaled  by  the 
reduced  x2  of  the  original  fit,  and  then  repeated  the  procedure. 
The  final  uncertainty  was  the  resulting  standard  deviation  of  the 
total  distribution. 

We  observed  HR  8799  with  multiple  calibrator  stars  every 
night  except  one  to  check  on  the  behavior  of  the  calibrators 
themselves.  For  example,  when  three  calibrator  stars  were  used, 
we  used  calibrator  1  to  determine  the  angular  diameters  of 
calibrator  2  and  calibrator  3  to  make  sure  the  stars  were  reliable. 
The  calibrator  HD  214698  shows  the  largest  scatter,  which  is 
expected  because  it  is  the  smallest  calibrator  of  the  three.  There 
will  naturally  be  more  uncertainty  when  measuring  its  calibrated 
visibilities  using  the  other  two  larger  stars.  The  scatter  seen  in 
the  night-to-night  measurements  is  expected  because  the  stars 
are  small  and  very  nearly  unresolved. 


3.2.  Stellar  Radius,  Luminosity,  and  Effective  Temperature 

At  a  distance  of  39.4  ±  1.1  pc  (van  Leeuwen  2007),  HR 
8799’s  0LD  of  0.342  ±  0.008  mas  corresponds  to  a  stellar 
radius9  of  1.44  ±  0.06  Rq,  corresponding  to  a  precision  of  4%. 
We  note  that  this  radius  is  8%  larger  than  the  radius  inferred 
from  photometric  and  temperature  considerations  in  Gray  & 
Kaye  (1999)  and  subsequently  adopted  in  recent  asteroseismic 
analyses  Moya  et  al.  (e.g.,  2010a,  2010b);  only  1%  of  this 
discrepancy  can  be  attributed  to  the  pre-revised  HIPPARCOS 
distance  used  in  the  calculation  of  Gray  &  Kaye  (1999). 

In  order  to  determine  the  luminosity  (L)  and  7'cff  of  HR  8799, 
we  first  constructed  its  photometric  energy  distribution  (PED) 
from  the  averages  (when  multiple  measurements  were  avail¬ 
able)  of  Johnson  UBV  magnitudes,  Stromgren  uvby  magni¬ 
tudes,  and  Two  Micron  All  Sky  Survey  (2MASS)  JHK  mag¬ 
nitudes,  as  published  in  Crawford  et  al.  (1966),  Breger  (1968), 
Eggen  (1968),  Mermilliod  (1986),  Schuster  &  Nissen  (1986), 
Mendoza  et  al.  ( 1 990),  Hauck  &  Mermilliod  (1998),  Gezari  et  al. 
(1999),  H0g  et  al.  (2000),  Cutri  et  al.  (2003),  and  Olsen  (1983, 
1993,  1994).  Table  5  summarizes  the  adopted  magnitudes;  the 
assigned  uncertainties  for  the  2MASS  infrared  measurements 


9  We  define  1  solar  radius  to  be  6.960  x  1010  cm,  consistent  with  an  average 

of  previous  measurements  (Emilio  et  al.  2012). 
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Figure  4.  #ld  fit  for  HR  8799  on  201 1  October  20.  The  symbols  are  the  same  as  in  Figure  1. 


are  as  reported,  and  for  optical  measurements  are  standard  devi¬ 
ations  of  multiple  measurements.  The  assigned  uncertainties  of 
the  optical  measurements  therefore  account  for  the  low  ampli¬ 
tude  optical  variability  (±0.02  mag)  observed  for  this  variable 
star  (Zerbi  et  al.  1999). 

The  bolometric  flux  (±bol)  of  HR  8799  was  determined  by 
finding  the  best-fit  (via  y2  minimization)  stellar  spectral  tem¬ 
plate  from  the  flux-calibrated  stellar  spectral  atlas  of  Pickles 
(1998).  This  best  PED  fit  allows  for  extinction,  using  the 
wavelength-dependent  reddening  relations  of  Cardelli  et  al. 
(1989).  The  best  fit  was  found  using  an  F0  V  template  with 
an  assigned  temperature  of  721 1  ±  90  K,  an  extinction  of  Ay  = 
0.00  ±  0.01  mag,  and  an  Tbol  of  1.043  ±  0.012  x 
10~7  erg  s'1  cm~2  (a  1.1%  precision).  Figure  1 1  shows  the  best 
fit  and  the  results  are  listed  in  Table  6.  To  check  for  possible 
systematic  biases  in  our  adopted  prescription,  we  also  fit  the 
PED  using  synthetic  Kurucz  spectral  models,10  assuming  no 
extinction.  The  resulting  Fbol  estimate  is  consistent  to  within 
the  1.1%  error  reported  above,  further  validating  our  technique 
and  measured  Fbol- 


10  Available  to  download  at  http://kumcz.cfa.harvard.edu 


The  bolometric  flux  was  then  combined  with  the  distance  to 
HR  8799  to  estimate  its  luminosity  (L  =  4tt  d2  /’boi.X  which 
yielded  a  value  of  5.05  ±  0.29  LQ.  The  uncertainty  in  the 
luminosity  (6%)  is  predominantly  set  by  the  uncertainty  in  the 
distance.  The  I'hoi.  was  also  combined  with  the  star’s  0\  n  to 
determine  its  effective  temperature  by  inverting  the  relation, 

±bol=  \etDaT^,  (2) 

where  a  is  the  Stefan-Boltzmann  constant.  This  produces  an 
effective  temperature  of  7203  ±  87  K,  determined  to  a  precision 
of  1%.  Because  /xld  is  chosen  using  a  given  7Q. n ,  we  used  our 
new  reff  value  to  select  /r ld  and  iterated,  /lld  increased  by  only 
0.01  to  0.49  ±  0.02,  0LD  increased  by  0.001  mas  to  0.342  ± 
0.008  mas,  and  reg-  decreased  by  10  K  to  7193  K.  The  very 
slight  change  in  0ld  did  not  affect  the  radius  calculation  at  all. 

We  note  that  this  Te g  is  nearly  identical  to  the  assigned 
temperature  of  the  best-fit  stellar  templates  used  to  calculate 
±bol-  However,  our  interferometric  Te ff  is  significantly  less 
model  dependent,  and  it  is  important  to  verify  the  accuracy 
of  the  Teff  predicted  by  the  PED.  Moreover,  our  measured  7’cn 
is  also  consistent  with  the  detailed  spectral-type  classification 
by  Gray  &  Kaye  (1999)  of  kA5  hFO  mA5  v  7.  Boo,  following 
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Figure  5.  #ld  fit  for  HR  8799  on  201 1  October  21.  The  symbols  are  the  same  as  in  Figure  1. 


the  spectral-type  temperature  scale  for  dwarf  stars  assembled  in 
Kraus  &  Hillenbrand  (2007).  The  hydrogen  lines  are  a  better 
tracer  of  the  stellar  Tt g-  (hFO;  7200  K);  the  metal  lines  suggest 
a  higher  Te g-  (mA5;  8200  K)  only  because  its  atmosphere  is 
slightly  metal  depleted. 

A  potential  bias  in  the  size  measurement  of  any  A-type  star 
is  oblateness  caused  by  rapid  rotation.  Royer  et  al.  (2007) 
measured  a  v  sin  /  of  49  km  s'  1  and  assuming  the  planets  orbit  in 
the  same  plane  as  stellar  rotation,  the  actual  rotation  velocity  will 
increase  from  49  to  104  km  s'  1  when  the  28°  inclination  is  taken 
into  account  (Soummer  et  al.  2011).  For  an  A5  V  star  with  an 
approximate  mass  of  2.1  MQ  (Cox  2000),  our  measured  radius 
of  1.44  Re,  and  the  relation  between  M,  oblateness,  and  v  sin/ 
described  in  van  Belle  et  al.  (2006),  the  predicted  oblateness  of 
HR  8799  is  only  ~2%  and  thus  within  the  lcr  errors  in  6*ld-  Using 
the  F0  V  spectral  type,  which  is  a  closer  fit  to  our  measured  Teg , 
the  estimated  mass  is  1.6  Me  and  the  predicted  oblateness  is 
2.6%,  still  a  small  effect. 

4.  DISCUSSION 

4.1.  Habitable  Real-estate  of  HR  8799 


have  not  yet  been  detected.  This  would  present  the  possibility  of 
life  if  the  planets  were  in  the  habitable  zone  (HZ)  of  the  star,  so 
we  used  our  new  precise  measurements  to  improve  the  estimate 
of  the  system’s  HZ.  We  used  the  following  equations  from  Jones 
et  al.  (2006): 

SbATeB)  =  (4.190  x  10“8  re2ff)  -  (2.139  x  10~4  Tcff)  +  1.296 

(3) 

and 

SbATeB )  =  (6.190  X  10~9  rjf)  -  (1.319  x  10"5  TeS ) 

+  0.2341,  (4) 

where  Sb,i(TeS )  and  Sb,0(TeB )  are  the  critical  fluxes  at  the  inner 
and  outer  boundaries  in  units  of  the  solar  constant.  The  inner 
and  outer  physical  boundaries  ri  o  in  AU  were  then  calculated 
using 

ri=n^  and  ra=f^.  (5) 

V  SbATetf)  V  Sb,o(Tes) 


HR  8799  is  currently  the  only  directly  imaged  multiple 
planet  system  and  there  may  be  smaller,  more  Earth-like  planets 
orbiting  the  star  or  even  moons  orbiting  the  imaged  planets  that 
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These  equations  assume  that  the  HZ  is  dependent  on  distance 
only  and  do  not  take  other  effects,  such  as  tidal  heating,  into 
account.  The  inner  boundary  is  the  limit  where  a  runaway 
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Figure  6.  $ld  fit  for  HR  8799  on  201 1  October  22.  The  symbols  are  the  same  as  in  Figure  1. 


greenhouse  effect  would  evaporate  any  surface  water  while  the 
outer  boundary  is  the  limit  where  a  cloudless  atmosphere  could 
maintain  a  surface  temperature  of  273  K.  Jones  et  al.  (2006)  note 
that  these  equations  produce  a  conservatively  small  HZ  and  the 
actual  HZ  may  be  wider. 

We  obtained  HZ  boundaries  of  1.62  AU  and  3.32  AU.  HR 
8799’s  planets  have  semimajor  axes  of  14.5-68  AU  (Marois  et  al. 
2008,  2010).  There  is  no  chance  the  planets  orbit  anywhere  near 
the  HZ  unless  they  are  highly  eccentric,  which  is  a  configuration 
more  likely  to  be  unstable. 

4.2.  Effective  Abundance  of  the  X  Boo  Stars  and  HR  8799 

The  accurately  determined  stellar  properties  of  HR  8799  also 
allow  for  detailed  comparisons  with  stellar  evolutionary  models 
from  which  a  mass  and  age  can  be  inferred.  However,  correctly 
interpreting  masses  and  especially  ages  depends  critically  upon 
knowing  the  internal  abundances  of  this  peculiar  abundance  star. 
As  noted  in  Section  1,  the  atmospheric  abundances  of  Fe-peak 
elements  are  distinctly  subsolar  ([Fe/H]  ~  —  0.4  dex)  while 
the  abundances  of  C,  N,  O,  and  S  are  essentially  solar.  Because 
these  surface  abundances  may  not  directly  trace  the  internal 
abundances,  it  is  unclear  which,  if  any,  of  the  available  uniformly 
scaled  abundance  models  to  adopt  for  these  comparisons; 


we  subsequently  refer  to  the  abundance  of  the  uniformly 
scaled  model  that  predicts  properties  most  consistent  with 
observational  constraints  as  the  effective  abundance. 

As  an  example  of  the  effect  on  the  inferred  ages,  compar¬ 
isons  of  HR  8799  with  models  that  assume  subsolar  abundances 
throughout  yield  ages  of  1.05  ±  0.26  Gyr  for  [Fe/H]=  —0.27 
or  1.71  ±  0.18  Gyr  for  [Fe/H]  =  -0.43,  while  those  that  as¬ 
sume  near-solar  abundances  ([Fe/H]  =  +0.05)  yield  an  age  of 
<0. 1  Gyr;  the  specific  physical  assumptions  used  in  these  mod¬ 
els  are  described  in  Section  4.3.  Obviously,  the  polyabundic 
atmospheres  of  X  Boo  stars  and  their  uncertain  internal  abun¬ 
dance  compromises  the  validity  of  these  comparisons  and  the 
inferred  values.  Nevertheless,  we  present  the  available  observa¬ 
tional  evidence  for  constraining  the  metallicity  of  X  Boo  stars, 
and  HR  8799  in  particular,  and  suggest  the  effective  abundance 
of  HR  8799,  and  perhaps  all  X  Boo  stars,  is  near  solar. 

The  observationally  favored  theory  to  explain  the  7,  Boo 
phenomenon,  originally  proposed  by  Venn  &  Lambert  (1990), 
is  that  the  atmospheres  have  been  polluted  by  the  accretion  of 
Fe-peak  depleted  gas.  Depletion  is  believed  to  occur  because 
of  grain  formation;  a  similar  depletion  of  Fe-like  elements  has 
been  observed  in  interstellar  clouds  (Morton  1974).  The  grains 
that  themselves  accrete  Fe-elements  are  consequently  inhibited 
from  accreting  on  to  the  star  because  of  the  stronger  radiation 
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Figure  7.  $ld  fit  for  HR  8799  on  201 1  September  30.  The  symbols  are  the  same  as  in  Figure  1. 


pressure  they  experience.  It  is  not  clear  why  this  accretion 
onto  the  grains  occurs,  or  where  the  accreting  material  comes 
from.  It  may  be  interstellar,  but  in  many  cases  appears  to  be 
circumstellar.  As  noted  by  Gray  &  Corbally  (2002),  all  four  X 
Boo  stars  within  40  pc  (including  HR  8799)  exhibit  an  infrared 
excess  that  is  interpreted  as  the  presence  of  circumstellar  dust, 
while  only  ~18%  of  noil-/.  Boo  A-type  stars  exhibit  such 
excesses.  The  apparent  depletion  of  these  elements  requires 
only  relatively  low  accretion  rates  (40  13  MQ  yr  1 ),  because 
of  the  shallow  convective  zones  in  A-type  stars  (Charbonneau 
1991).  Accretion  at  these  rates  will  quickly  establish  abundance 
anomalies  within  a  few  Myr,  but  these  anomalies  will  likewise 
disappear  in  as  little  as  1  Myr  once  the  accretion  has  terminated 
(Turcotte  &  Charbonneau  1993);  the  interesting  implication  is 
that  all  X  Boo  stars  are  either  currently  accreting,  or  have  only 
recently  terminated  their  accretion.  If  this  favored  theory  is 
correct,  it  suggests  that  the  depleted  Fe-peak  abundances  do 
not  represent  the  internal  abundances  of  these  stars.  Given  this 
and  the  typical  limiting  main-sequence  lifetime  of  ~2  Gyr  for 
these  intermediate  mass  Population  I  stars,1 1  one  would  expect 
their  internal  abundance  to  be  on  average  close  to  solar. 

We  note  that  an  alternative  mechanism  to  explain  the  X  Boo 
phenomenon  proposed  by  Michaud  &  Charland  (1986)  suggests 

1 1  Studies  have  shown  that  stars  currently  forming  in  open  clusters  are 
generally  metal-rich  (e.g.,  Santos  et  al.  2009). 


that  the  depletion  of  Fe-peak  elements  is  a  result  of  diffusion 
and  mass  loss.  In  this  case,  X  Boo  stars  are  much  closer  to 
the  end  of  their  main-sequence  lifetimes.  To  be  effective,  the 
star  would  need  to  be  losing  mass  at  a  rate  of  the  order  of 
I  (A  1 3  M0  yr~'  forafew  times  108  yr  and  thus  implies  that  the! 
Boo  phenomenon  be  restricted  to  the  end  of  the  main-sequence 
evolutionary  phase  for  A  stars.  This  predicted  timescale  is 
difficult  to  reconcile  with  the  discovery  of  very  young  X  Boo 
stars  in  the  Orion  OBI  associations  (age  <10  Myr),  and  the 
lack  of  any  X  Boo  stars  in  intermediate  age  open  clusters  (Gray 
&  Corbally  2002).  Evidence  such  as  the  higher  fraction  of 
these  stars  with  circumstellar  debris  disks,  noted  above,  instead 
suggest  that  the  X  Boo  phenomenon  is  more  common  among 
young  A  stars.  If  true,  this  would  strengthen  the  case  that 
the  internal  abundances  of  these  Population  I  stars,  at  least  on 
average,  should  be  close  to  solar.  In  particular,  we  note  that  the 
existence  of  X  Boo  stars  in  the  Orion  star-forming  region  implies 
that  their  primordial  abundances  are  solar,  consistent  with  other 
stars  in  this  cluster  (Nieva  &  Przybilla  2012).  If  the  depletion 
of  Fe-like  elements  is  restricted  to  only  the  outer  atmospheres, 
then  one  can  conclude  that  in  these  cases  the  internal  abundances 
should  likewise  be  solar. 

As  described  in  Marois  et  al.  (2008),  the  space  motions  of 
HR  8799  are  consistent  with  those  of  young  stars  in  the  solar 
neighborhood.  We  investigated  this  further  by  comparing  the 
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Figure  8.  #ld  fit  for  HR  8799  by  calibrator.  The  symbols  are  the  same  as  in  Figure  1 . 


UVW  space  motions  of  HR  8799  to  those  of  nearby  moving 
groups  assembled  in  Torres  et  al.  (2008);  kinematic  association 
with  a  moving  group  would  not  only  help  constrain  the  age  of 
the  system  but  also  the  primordial  abundance  of  the  star.  We 
calculated  UVW  space  motions  using  HR  8799’s  HIPPARCOS 
distance  and  parallax  (van  Leeuwen  2007),  and  a  radial  velocity 
of  — 12.6  ±  1.3  km  s-1  from  Barbier-Brossat  &  Figon  (2000). 
This  yielded  velocities  of  U  =  —12.24  ±  0.37  km  s_1,  V  — 
-21.22  ±  1.10  km  s_1,  and  W  =  -7.15  ±  0.86  km  s"1;  these 
are  all  within  ~1  km  s-1  of  the  values  reported  in  Marois 
et  al.  (2008).  Figure  12  illustrates  the  space  motion  of  HR  8799 
relative  to  that  of  several  kinematically  similar  moving  groups 
after  adopting  moving  group  velocities  and  uncertainties  from 
Torres  et  al.  (2008).  Although  the  space  motion  of  HR  8799  is 
not  consistent  to  within  1  a  of  any  of  these  groups,  it  is  consistent 
to  within  3er  of  two  of  these  groups:  Columba  (1.2er),  and  e  Cha 
(2.2(7).  We  thus  confirm  the  report  by  Marois  et  al.  (2008)  of 
the  potentially  youthful  kinematics  of  HR  8799,  which  favors 
but  cannot  confirm  an  age  less  than  <1  Gyr,  and  likely  a  solar 
abundance  consistent  with  many  of  these  stars. 

As  described  in  Section  1,  there  are  two  other  aspects  of  HR 
8799  that  argue  for  a  young  age  for  the  star  and  system.  The 
strong  infrared  excess  of  HR  8799  statistically  favors  an  age 


less  than  ~500  Myr  (Gaspar  et  al.  2009),  despite  the  challenge 
raised  by  Moya  et  al.  (2010a).  An  even  stronger  youthful  age 
restriction  is  imposed  by  the  dynamical  considerations  of  Moro- 
Martfn  et  al.  (2010).  Although  their  conclusions  require  use 
of  the  predictions  of  poorly  constrained  planetary  evolutionary 
models,  even  considerable  uncertainties  in  the  mass  estimates 
require  the  age  of  the  system  to  be  less  than  a  few  100  Myr, 
otherwise  the  planets  would  be  too  massive  to  be  orbitally 
stable.  The  ages  much  less  than  1  Gyr  for  HR  8799  can  only  be 
reconciled  if  the  adopted  effective  abundances  are  near  solar. 

Recently,  Moya  et  al.  (2010b)  used  the  y  Dor-type  pulsations 
of  HR  8799  to  asteroseismically  constrain  the  star’s  internal 
metallicity.  They  find  a  best-fit  internal  metallicity  of  —0.32  ± 
0.1,  assuming  that  the  rotation  axis  of  the  star,  which  affects 
the  rotation  speed  of  the  star,  is  somewhat  highly  inclined 
relative  to  our  line  of  sight  (50°).  This  high  inclination  is 
inconsistent  with  the  more  face-on  orientation  of  HR  8799’s  disk 
(Su  et  al.  2009)  and  the  apparent  orbital  plane  of  its  companions 
(Lafreniere  et  al.  2009),  although  these  need  not  be  coplanar. 
The  asteroseismology  analysis  is  compromised  by  the  use  of  an 
assumed  radius  that  is  too  small  by  8%  and  a  temperature  that 
is  too  hot  by  3%  relative  to  the  precisely  determined  values 
in  this  study.  The  measured  R  and  Tctj  presented  here  will 
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Figure  9.  #ld  fit  for  HR  8799  using  all  calibrators.  The  symbols  are  the  same  as  in  Figure  1. 


be  valuable  to  constrain  the  model  parameter  space  in  future 
asteroseismic  modeling  efforts,  which  will  possibly  constrain 
the  internal  metallicity. 

Moya  et  al.  (2010b)  do  find  possible — though  less 
likely — solutions  with  internal  abundances  close  to  solar 
([Fe/H]  =  —0.12)  for  closer  to  pole-on  orientations;  this  is 
the  most  metal-rich  metallicity  reported  in  their  study.  Overall, 
we  find  the  results  of  this  metallicity  analysis  to  be  inconclusive 
primarily  because  of  the  large  uncertainty  in  the  inclination  of 
the  star’s  rotation  axis,  which  Moya  et  al.  acknowledge  is  the 
limiting  factor  in  their  analysis.  If  anything,  the  result  that  they 
find  some  acceptable  solutions  with  near-solar  metallicities  and 
young  ages  corroborated  previous  asteroseismology  studies  of 
X  Boo  stars.  Using  densities  inferred  from  stellar  pulsations, 
Paunzen  et  al.  (1998)  found  the  location  of  most  X  Boo  stars 
showing  <5  Scuti  pulsations  on  a  plot  of  average  density  versus 
period  to  be  consistent  with  the  positions  of  8  Scuti  stars  with 
solar  abundances.  They  interpret  this  as  evidence  that  the  low  Fe- 
peak  abundances  are  restricted  to  the  surface  of  X  Boo  stars  and 
do  not  represent  the  state  of  the  interior  composition.  Although 
the  variations  in  HR  8799  are  compatible  with  y  Dor  pulsations, 
the  large  number  of  discovered  hybrid  8  Scu-y  Dor  pulsators 
(e.g.,  Rowe  et  al.  2006;  Grigahcene  et  al.  2010)  suggests  that 
this  argument  could  be  applicable  for  all  X  Boo  stars. 

In  summary,  based  on  the  observational  evidence  supporting 
the  accretion  of  clean  gas  theory  explaining  the  abundance 


pattern  of  X  Boo  stars,  the  existence  of  at  least  some  X 
Boo  stars  in  an  OB  Association,  the  evidence  of  youth  for 
many  X  Boo  stars  (and  especially  HR  8799),  the  restriction 
of  ages  less  than  ~2  Gyr  for  Population  I  stars,  and  the 
available  asteroseismic  constraints  on  internal  metallicity,  we 
conclude  that  the  most  appropriate  effective  abundances  for  HR 
8799,  and  quite  possibly  all  X  Boo  stars,  is  near  solar.  Until 
more  sophisticated  evolutionary  models  that  account  for  the 
polyabundic  atmospheres  and  possibly  interiors  of  these  stars 
are  developed,  we  are  hopeful  that  X  Boo  stars  in  clusters  or  with 
lower  mass  companions  will  be  discovered,  which  will  enable 
improvements  in  both  assigning  effective  abundances  and  tests 
of  the  validity  of  using  uniformly  scaled  abundances. 

4.3.  Mass  and  Age  of  HR  8799  and  Implications 
for  Its  Companions 

Following  the  above  arguments,  we  estimated  the  mass  and 
age  of  HR  8799  using  the  stellar  evolutionary  models  of  the 
Yonsei-Yale  group  (Y2;  Yi  et  al.  2001),  updated  by  Demarque 
et  al.  (2004)  to  account  for  convective  core  overshoot.  The 
models  handle  convection  using  mixing  length  prescription, 
adopting  a  mixing  length  of  1.7432  times  the  pressure  scale 
height  that  is  set  by  comparisons  with  a  solar  model.  The  model 
uses  the  solar  abundances  of  Grevesse  et  al.  (1996)  and  OPAL 
radiative  opacities  (Rogers  &  Iglesias  1995;  Iglesias  &  Rogers 
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Figure  10.  Comparison  of  #ld  fits  by  individual  night  and  calibrator.  The  vertical  solid  line  represents  the  final  adopted  #ld  (0.342  ±  0.008  mas,  the  top  line)  and  the 
vertical  dashed  lines  are  the  errors  in  that  fit.  Table  4  lists  the  numerical  values.  Note  the  reduced  scatter  in  the  201 1  data  vs.  the  2010  data. 
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Figure  11.  HR  8799  PED  fit.  Upper  panel:  the  solid-line  spectrum  is  an  F0  V  spectral  template  from  Pickles  (1998).  The  crosses  indicate  photometry  values  from  the 
literature.  The  horizontal  bars  represent  bandwidths  of  the  filters  used.  The  X-shaped  symbols  show  the  flux  value  of  the  spectral  template  integrated  over  the  filter 
transmission.  Lower  panel:  the  crosses  are  the  residuals  around  the  fit  in  fractional  flux  units  of  photometric  uncertainty.  For  more  details,  see  Section  4. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


1996)  for  the  interior  and  the  Alexander  &  Ferguson  (1994) 
opacities  for  the  cooler,  outer  regions  of  the  star.  Additional 
information  regarding  the  input  physics  assumed  in  these  models 
can  be  found  in  Yi  et  al.  (2001)  and  Demarque  et  al.  (2004). 


However,  since  a  set  solar  metallicity  models  is  not  provided, 
we  followed  the  recommendation  of  the  modelers  themselves 
and  generated  a  set  of  solar  metallicity  models  by  linearly 
interpolating  isochrones  and  mass  tracks  between  nearest  the 
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Figure  12.  UVW  space  velocities  of  HR  8799  (black  point)  with  1  a  (error  bars)  and  2a  (gray  ellipse)  errors.  Also  plotted  are  the  UVW  space  velocities  and  lcr  errors 
(as  ellipses)  of  four  young  stellar  associations  with  similar  kinematics,  as  taken  from  Torres  et  al.  (2008).  HR  8799  matches  the  space  velocity  of  the  roughly  30  Myr 
old  Columba  association  to  1 ,2cr . 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 
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Figure  13.  Radius  and  temperature  of  HR  8799  plotted  along  with  isochrones 
(blue  lines)  and  mass  tracks  (solid  red  lines)  from  Demarque  et  al.  (2004),  with 
solar  abundances.  Note  that  the  mass  tracks  predict  that  these  stars  are  still 
gravitationally  settling  at  0.02  Gyr,  reach  their  smallest  size  at  0.04  Gyr,  and 
expand  in  size  thereafter. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


metallicity  models  of  [Fe/H]  =  —0.27  dex  and  [Fe/H]  =  +0.05 
dex  (models  x74z01  and  x71z04  models,  respectively).  The  solar 
metallicity  models  are  illustrated  in  Figure  13. 

At  HR  8799’s  temperature,  stars  are  predicted  to  have  fully 
contracted  to  their  smallest  size,  or  the  zero-age  main  sequence 
(ZAMS)  at  40  Myr,  and  then  begin  expanding  due  to  stellar 
evolutionary  effects  at  a  much  slower  rate  afterward.  As  such, 
HR  8799  is  either  contracting  onto  the  ZAMS  or  expanding  from 
it,  and  we  use  the  40  Myr  isochrone  as  the  upper  limit  for  the 
age  in  the  pre-main-sequence  scenario  and  as  the  lower  limit  for 
the  age  in  the  post-main-sequence  scenario.  More  specifically, 
we  find  that  if  HR  8799  is  contracting  onto  the  ZAMS,  it  has 
an  age  of  33/.713i  Myr  and  a  mass  of  1.516^^  Mq.  If  HR 
8799  is  expanding  from  the  ZAMS,  we  find  it  to  have  an  age 
of  90*3581  Myr  and  a  mass  of  1.51 3/qq2234  Mq.  These  masses  and 
ages  are  found  by  interpolating  between  the  solar  metallicity 
models  described  above.  The  masses  are  consistent  with  the 
mass  used  in  Marois  et  al.  (2008),  which  was  1.5  ±  0.3  Mq. 
We  remind  the  reader  that  these  quoted  errors  on  the  star’s  mass 


and  age  are  statistical  and  therefore  do  not  take  into  account 
uncertainties  in  the  metallicity  and/or  models  themselves. 

With  an  age  of  <0.1  Gyr,  the  companions  that  HR  8799 
harbors  are  more  likely  to  have  planetary  masses.  As  explained 
in  Marois  et  al.  (2008),  planetary  mass  objects  could  only 
be  as  bright  as  these  observed  companions  are  if  they  are 
young.  That  inferred  age,  and  thus  companion  masses,  depends 
critically  upon  adopted  evolutionary  model  abundance.  While 
near  solar  seems  likely,  even  slightly  subsolar  abundances  can 
give  uncertainties  in  age  that  extend  to  above  ~500  Myr, 
increasing  the  inferred  companion  masses  by  at  least  a  factor  of 
two. 

5.  SUMMARY 

We  measured  the  angular  diameter  of  HR  8799  using  the 
CHARA  Array  interferometer  and  used  our  new  value  of 
0.342  ±  0.008  mas  to  calculate  the  star’s  physical  radius 
(1.44  ±  0.06  Rq),  luminosity  (5.05  ±  0.29  L0),  and  effective 
temperature  (7 1 93  ±  87  K)  by  combining  our  measurement  with 
information  from  the  literature.  We  used  our  7/n  measurement 
to  determine  the  size  of  the  HZ,  which  is  well  inside  the  orbits 
of  any  of  the  companions  detected  to  date. 

Based  on  a  variety  of  techniques,  we  concluded  that  the  most 
appropriate  abundances  for  HR  8799  are  close  to  solar.  We 
combined  our  R  and  Teg  values  with  Y2  isochrones  to  estimate 
the  star’s  mass  and  age  in  two  scenarios:  1.516^^  Mq  and 
33/7i3  2  Myr  if  the  star  is  contracting  onto  the  ZAMS  or 
1.513^00^4  Mo  and  90*35801  Myr  if  it  is  expanding  from  it.  In 
either  case,  this  young  age  implies  the  imaged  companions 
are  planets  and  not  brown  dwarfs.  The  only  case  in  which  the 
companions  would  be  close  to  brown  dwarf  mass  is  if  the  highest 
age  of  471  Myr  was  the  true  case,  and  even  then  the  masses 
would  be  on  the  exoplanet/brown  dwarf  cusp. 

We  thank  Gerard  van  Belle  for  his  insight  on  the  nature 
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